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1. Introduction 


Nanotechnology stands for a new paradigm in materials science that aims at exploring 
and controlling new physical phenomena and enhanced functionalities that emerge at the 
nanoscale. These opportunities are a direct consequence of reduced dimensionality and/or 
interfacial phenomena from proximity effects between dissimilar materials Zubko et al. (2011). 
Progress in growth techniques Chambers (2010); Eckstein & Bozovic (1995); Martin et al. 
(2010); McKee et al. (1998); Posadas et al. (2007); Reiner et al. (2009); Schlom et al. (1992); 
Vaz et al. (2009a); Vrejoiu et al. (2008), nanoscale characterization tools Zhu (2005), and first 
principles calculations Cohen (2000); Fennie (2008); Picozzi & Ederer (2009); Rabe & Ghosez 
(2007); Spaldin & Pickett (2003); Waghmare & Rabe (2005) have been instrumental to our 
present ability to control matter down to the atomic scale and to fabricate nanoscale device 
structures with the potential for technological applications. Examples of current research 
work that aims at addressing some of the current grand challenges include the search for 
ultrasensitive sensors and actuators for applications in areas such as medicine and energy 
harvesting, the development of smaller and more energy efficient electronic devices that 
could replace current CMOS switches, and the design of intelligent systems that incorporate 
complex operations at the core processing level. While the approach employed to date has 
relied on building up complexity from basic building blocks (e.g., complex microprocessor 
units formed of MOSFET devices), a new approach is being developed that directly explores 
the complex state of matter to achieve integrated functionalities and more complex operations 
at a fundamental level. Key to this effort has been the sustained research work aimed 
at understanding the properties of strongly correlated systems, and at controlling such 
properties down to the atomic scale, from a device physics perspective. 

In this context, a particularly interesting class of materials are so-called multiferroic systems, 
which are characterized by the presence of simultaneous magnetic and ferroelectric order. As 
such, they offer the possibility of achieving control of the magnetic state via applied electric 
fields, or vice versa, which could find applications in ultrasensitive magnetic sensors and 
transducers, solid state transformers, magnetoelectrooptic devices, energy harvesting and 
storage, and new spin-based logic devices in the context of spintronics Bibes & Barthélémy 
(2007); Cibert et al. (2005); Žutić et al. (2004). In most single phase multiferroic materials, the 
origin of magnetic and ferroelectric orders is largely independent, with the consequence that 
the coupling between magnetism and ferroelectricity (mediated by the spin-orbit coupling) 
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is weak; in those instances where the coupling is strong, the critical temperatures tend to be 
small Khomskii (2006; 2009); Picozzi & Ederer (2009). To overcome the weak magnetoelectric 
coupling of single-phase multiferroic materials, alternate approaches have been developed 
that explore proximity and interfacial effects between magnetic and ferroelectric materials to 
form composite structures with enhanced coupling between electric and magnetic properties. 
By judiciously engineering the interfacial properties at the nanoscale, a strong coupling 
between magnetic and ferroelectric order parameters can be achieved. This new class of 
artificially structured composite materials exhibits magnetoelectric couplings that are orders 
of magnitude larger that those typical of single-phase, intrinsic multiferroics Fiebig (2005); Ma 
et al. (2011); Vaz et al. (2010a). 

An example of the promise afforded by this approach is provided by the particular case 
of the multiferroic perovskite BiFeO3, characterized by magnetic and ferroelectric critical 
temperatures well above room temperature (T® = 643 K and TF = 1100 K, respectively) 
Catalan & Scott (2009). BiFeO3 has generated much interest recently, following the first report 
of the growth of epitaxial thin films Wang et al. (2003) and the demonstration of very large 
electric polarizations in high quality single crystalline films and in bulk crystals Lebeugle et al. 
(2007); Shvartsman et al. (2007); Wang et al. (2003). In this system, ferroelectricity originates 
from the lone-pair active Bi cations Neaton et al. (2005); Ravindran et al. (2006), while the 
magnetic order originates from the oxygen-mediated superexchange interaction between 
the Fe cations, which favors an antiferromagnetic coupling between nearest neighbor spins 
Kiselev et al. (1963). The magnetic state is modified further by a break in center of symmetry 
and the presence of a ferroelectric polarization, which gives rise to a local spin canting between 
the two spin sublattices (and to a weak magnetic moment) through the Dzyaloshinskii-Moriya 
interaction Dzialoshinshkii (1957); Moriya (1960). In addition, the coupling of the polarization 
to gradients of the magnetization leads to an inhomogeneous spin configuration characterized 
by an incommensurate rotation of the total local spin along a (101)p¢ direction (indexed 
to the pseudocubic perovskite structure) and lying in the {121}pc plane, defined by the 
cycloid propagation direction and the electric polarization (with easy axes along (111)pc), 
with a period of about 62 nm (spin cycloid) Catalan & Scott (2009); Picozzi & Ederer (2009); 
Sosnowska et al. (1982). This spin cycloid averages out the magnetic moment and leads to a 
vanishing linear magnetoelectric coupling and to a small effective magnetoelectric response; 
however, at the nanoscale, there is a strong coupling between the electric polarization and the 
magnetic spins, since they are constrained to point perpendicular to each other, indicating 
that a change in the orientation of the electric polarization will result in a change in the 
spin direction Cazayous et al. (2008); Lebeugle et al. (2008). Such a phenomenon has been 
demonstrated experimentally Lee et al. (2008); Zheng et al. (2006) and has been explored in 
exchange-bias coupled multiferroic heterostructures to change the magnetization direction of 
a ferromagnetic layer exchange-coupled to the BiFeO3 Chu et al. (2008); Wu et al. (2010). This 
approach to magnetoelectric coupling illustrates the current trend towards engineering larger 
magnetoelectric couplings by relying on interfacial effects between different materials, an 
approach that can be traced back to the 1970s, when the first attempts to grow strain-mediated 
ferroelectric-ferromagnetic composites were made van Suchtelen (1972). In fact, the most 
common approach to date for achieving a magnetoelectric coupling in composites relies on 
strain to mediate the magnetic and electrical properties by inducing crystal deformations on 
either the magnetic phase through the converse piezoelectric effect, or in the ferroelectric 
phase through magnetostriction Fiebig (2005); Ma et al. (2011); Ramesh & Spaldin (2007); 
Thiele et al. (2007); Vaz et al. (2009b; 2010a). The effect is indirect, but can be optimized to yield 
large magnetoelectric responses by a suitable choice of the material components and device 
geometry Nan et al. (2008); Srinivasan (2010); Wang et al. (2009). Another promising route 
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involves a direct, charge-mediated magnetoelectric coupling in ferroelectric/ferromagnetic 
oxide composite multiferroic heterostructures, where the spin state of the magnetic oxide 
is controlled via charge doping induced by the electric polarization of a ferroelectric 
Molegraaf et al. (2009); Vaz et al. (2010b). For Pb(Zro.>Tio.s)O3 / Lag. g5r9.2MnO3 (PZT /LSMO) 
heterostructures, the effect is electronic in origin and results from a change in the magnetic 
spin configuration as a function of the ferroelectric polarization direction, demonstrating 
electric field control of magnetism in this system (Section 3). 

In this chapter, we consider the recent developments in the electric field control of magnetism 
in artificial heterostructures based on electrostatic doping. This approach has been explored in 
various systems, including dilute magnetic semiconductors, transition metals, and complex 
oxides; an overview of the work carried out in each of these systems is given in Section 2. 
The focus will be on complex oxide heterostructures, and in particular on our recent work 
demonstrating a strong magnetoelectric coupling in PZT/LSMO multiferroic heterostructures 
(Section 3). Due to strong electron correlations, complex oxides offer an inexhaustible range 
of possibilities for the study of novel phenomena that arise from the sensitivity of these 
materials to charge, strain, electric and magnetic fields, among other control parameters, with 
the attendant promise for device applications Imada et al. (1998); Tokura (2006); Tokura & 
Nagaosa (2000). 


2. Electrostatic control of magnetism in artificial heterostructures 


One approach to artificial multiferroic structures exploits the electric field effect to achieve 
an electrostatic modification of the charge carrier density and to induce changes in the 
magnetic state. The working concept is similar to that of a field effect transistor, where an 
induced or spontaneous electric polarization at the gate dielectric interface is screened by 
charge carriers from the channel layer, leading to charge accumulation or depletion over a 
thickness determined by the screening length of the material. In materials where the magnetic 
properties are intimately linked to charge, a change in carrier doping results in a change in 
the magnetic properties. The amount of charge carrier modulation required will depend on 
the particular system. For strongly correlated oxides, where typical carrier densities are of 
the order of 10°! cm~%, the requisite modulation in the charge carrier doping can be achieved 
by using a ferroelectric for the gate dielectric, an approach termed the ferroelectric field effect 
Ahn et al. (2006); Venkatesan et al. (2007). In this approach to electrostatic doping, the charge 
carriers screen the large surface bound charge of the ferroelectric; for a ferroelectric such as 
Pb(Zr,Ti)O3 (PZT), the charge carrier modulation is of the order of 1014 cm~?, much larger than 
is possible to attain using silicon oxide as the gate dielectric Ahn et al. (2003). The field effect 
approach has been explored to control a variety of properties in complex systems, including 
superconductivity Ahn et al. (1999); Caviglia et al. (2008); Frey et al. (1995); Parendo et al. 
(2005); R. E. Glover & Sherrill (1960); Talyansky et al. (1996) and metal-insulator transitions 
Dhoot et al. (2009); Hong et al. (2005; 2003). Ferroelectric gates have also been proposed 
for non-volatile field effect transistors (FET), where the on/off states are maintained by the 
ferroelectric polarization Brown (1957); Looney (1957); Miller & McWhorter (1992); Park et al. 
(2003). Electrostatic control of magnetism has been reported in the last decade for a variety 
of systems, including dilute magnetic semiconductors (DMS), transition metal ferromagnets, 
and complex oxides, as discussed briefly below. 

The phenomenology of the magnetoelectric coupling in materials has been discussed 
extensively, and we refer to a recent review for details and for the relevant literature Vaz et al. 
(2010a). The figure of merit that characterizes the coupling between the electric and magnetic 
order parameters is the magnetoelectric susceptibility, which measures the change in magnetic 
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moment for a given applied electric field. In the case where the magnetoelectric response is 
linear, the magnetoelectric susceptibility a = uodM/dEo (in S.I., where pg is the permeability 
of vacuum, M is the magnetization, and Eg is the external applied electric field) is well defined, 
but it is less so in the more general case where the magnetoelectric response is non-linear. In 
such instances, it is common to define an effective magnetoelectric constant corresponding 
to the change in magnetization for a given applied electric field (or conversely, the change 
in electric polarization for a given applied magnetic field). In composite systems relying 
on interfacial effects, it is useful to define a surface (interface) magnetoelectric coefficient &s, 
corresponding to the change in surface magnetization for a given applied electric field Duan 
et al. (2008); Fechner et al. (2008); Niranjan et al. (2009). With the definition given above, 
the linear magnetoelectric constant has units of s m! in S.I. units, while in cgs units it is 
dimensionless Rivera (1994); they are related one to the other by the speed of light in vacuum, 
such that a dimensionless relative magnetoelectric constant («,) independent of the system of 
units can be defined, in analogy with the magnetic and electric relative permitivities Hehl et al. 
(2008; 2009). Often, however, « is given in mixed units, such as Oe cm V-t, We list in Table 1 
the magnetoelectric response of charge-mediated multiferroic heterostructures reported in the 
literature, both in terms of the interfacial and relative magnetoelectric coupling coefficients. 


System 10a 101%; a, T(K) Reference 

bcc Fe(001) 0.002 0.024 0.0005 Theory Duan et al. (2008) 
hep Co(0001) 0.0008 0.016 0.0002 Theory Duan et al. (2008) 
fee Ni(001) 0.002 0.03 0.0005 Theory Duan et al. (2008) 

9 ML Fe/MgO 0.008 0.11 0.0023 Theory Niranjan et al. (2010) 


2 ML Fe/BaTiO3 10 200 3.0 Theory Duan et al. (2006) 
1 ML Fe/BaTiO3 16 230 48 Theory Fechner et al. (2008) 


Fe/BaTiO3/Pt 3 43 0.9 Theory Cai et al. (2009) 
1 ML Fe/PbTiO3 73 1000 22 Theory Fechner et al. (2008) 
Ni/BaTiO3/Pt 15 260 4.5 Theory Cai et al. (2009) 


hep Co/BaTiO3/Pt 4 81 1.2 Theory Cai et al. (2009) 

CrO, /BaTiO3/Pt 10 150 3.0 Theory Cai et al. (2009) 
Fe304/BaTiO3 20 200 57 Theory Niranjan et al. (2008) 
SrRuO3/SrTiO3 0.05 2 0.015 Theory Rondinelli et al. (2008) 
SrRuO3/BaTiO3 1.1 42 0.32 Theory Rondinelli et al. (2008) 
SrRuO3/BaTiO3 5.9 230 1.8 Theory Niranjan et al. (2009) 
PZT/LSMO (x = 0.2) 0.8 310 24 100K Molegraaf et al. (2009) 
PZT/LSMO (x = 0.2) 6.2 2900 22 100K Vaz et al. (2010b;c) 
PZT/LSMO (x = 0.2) -13.5 -6300 -49 180K Vaz etal. (2010c) 


Table 1. Values of the magnetoelectric coupling coefficient reported in the literature for 
charge-driven multiferroic heterostructures. « is given in units of Oe cm VTL, ws in units of 
Oe cm? V~!; T is the temperature. When not directly provided, the surface magnetoelectric 
coupling coefficient «s is estimated by multiplying « by 1 ML of the corresponding magnetic 
material; for the case of the PZT/LSMO structures, a; and a; are estimated by multiplying 
the experimental value of a by the LSMO film thickness. 


The nature of the magnetoelectric effect due to charge screening can be distinguished between 
(i) enhanced spin imbalance at the Fermi level due to screening and the corresponding 
modification in the magnetic moment of the system as a function of the electric field Zhang 
(1999); (ii) changes in magnetic moment due to changes in electronic bonding at the polarized 
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dielectric interface Duan et al. (2006); (iii) changes in the magnetic order with the charge 
density Gerhard et al. (2010); Kudasov & Korshunov (2007); Ovchinnikov & Wang (2008); 
Sun et al. (2010); Vaz et al. (2010b), whereby the magnetic state of the system is modified due 
to changes in the charge carrier density, either between magnetic and non-magnetic states, or 
between states with different magnetic spin configurations; and (iv) changes in the magnetic 
anisotropy that lead to different global magnetic states for different applied electric fields 
Maruyama et al. (2009); Niranjan et al. (2010). 


2.1 Electrostatic control of magnetism in dilute magnetic semiconductors 

The field effect approach to controlling magnetism is well suited for dilute magnetic 
semiconductors (DMS), such as (In,Mn)As, (Ga,Mn)As, and Mn;Ge;_,, where integration 
with semiconductor substrates, such as GaAs, follows naturally from the film growth process. 
In the DMS systems, the ferromagnetic interaction between the Mn spins (mediated by 
hole carriers) competes with the antiferromagnetic superexchange, and becomes dominant 
at sufficiently high hole doping Dietl et al. (2000). The demonstration of electric field 
modulation of magnetism in a 5 nm epitaxial (In,Mn)As layer was reported by Ohno et al. 
(2000), using a thick polyimide layer as the gate dielectric. By varying the charge carrier 
doping through the application of a gate voltage, a change in the critical temperature by 
about 2 K (Te = 25 K at zero electric field) is achieved between the accumulation and 
depletion states. The size of the effect is found to agree with a Zener model used to describe 
the onset of magnetism in zinc-bled magnetic semiconductors Dietl et al. (2000). Besides 
(In,Mn)As Chiba et al. (2003); Ohno et al. (2000), electric field control of magnetism in magnetic 
semiconductors has been reported for Mn,Geı_, Chen et al. (2007); Park et al. (2002); Xiu 
et al. (2010), (Zn,Mn)Se Kneip et al. (2006), and (Ga,Mn)As Chiba et al. (2008); Endo et al. 
(2010a); Nazmul et al. (2004); Owen et al. (2009); Riester et al. (2009); Stolichnov et al. (2008). 
In the report by Stolichnov et al. (2008), a ferroelectric polymer (polyvinylidene fluoride with 
trifluoroethylene, or P(VDF-TrFE)) is employed as the gate dielectric to achieve non-volatile 
control of ferromagnetism in a (Ga,Mn)As channel layer, as manifested by changes in the 
coercivity of the magnetic hysteresis loop and in the critical temperature as a function of the 
ferroelectric polarization direction. 


2.2 Electrostatic control of magnetism in transition metals 

The electrostatic control of magnetism in transition metals has been demonstrated by Weisheit 
et al. (2007), who report a modulation of the magnetic coercivity in ordered FePt and FePd 
intermetallic alloys subject to an applied field when immersed in an electrolyte. More 
recently, Maruyama et al. (2009) have shown that the perpendicular magnetic anisotropy of 2-4 
monolayers (ML) Fe films can be modified by electric fields by using a polyimide layer as the 
gate dielectric (up to 40% for applied fields of the order of 1 MV/cm). The effect is particularly 
promising and could be explored to devise Fe/MgO/Fe magnetic tunnel junctions, where 
the tunnel magnetoresistance (TMR) is tuned electrostatically by controlling the magnetic 
anisotropy—and the relative magnetization alignment—of the Fe layers. Record high tunnel 
magnetoresistance ratios have been reported for Fe/MgO/Fe Bowen et al. (2001); Parkin et al. 
(2004); Yuasa et al. (2004), making this system a promising candidate for spintronic device 
applications. Indeed, control of the magnetic easy axis from in-plane to out of plane has been 
recently demonstrated for FeggCo29 Shiota et al. (2009), while modulation of the magnetic 
anisotropy has been shown in Feg9Co29 /MgO/Fe tunnel junctions Nozaki et al. (2010) and in 
MgO/CoaoFe40B20/ Ta structures Endo et al. (2010b). Another demonstration of electric field 
modulation of magnetism relies on the magnetic and structural instabilities of elemental Fe, 
where the application of an electric field is found to induce changes in the crystal structure of 
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Fe/Cu(111) islands, from the ferromagnetic bcc to the antiferromagnetic fcc structure Gerhard 
et al. (2010); first principles calculations show that the effect originates from changes in the Fe 
interatomic distances with the applied electric field as a result of electronic charge screening 
that tilts the energy of the system to favor a ferromagnetic or antiferromagnetic state Gerhard 
et al. (2010). 

The theoretical study of the magnetoelectric coupling in metal-based multiferroic 
heterostructures has been the subject of intensive investigation. One can distinguish two 
approaches, one that considers the effect of charge screening in free standing metal layers 
subject to an external electric field, and the other that considers charge screening at the 
interface with a dielectric. In the latter case, contributions from chemical bonding between 
metal and dielectric need to be considered. The nature of the magnetoelectric effect that 
appears in metals can also be distinguished between enhanced spin imbalance at the Fermi 
level due to screening Zhang (1999), changes in the chemical bonding at the interface with the 
polarized dielectric, changes in the magnetic order with the charge density Gerhard et al. 
(2010); Kudasov & Korshunov (2007); Ovchinnikov & Wang (2008); Sun et al. (2010), and 
changes in the magnetic anisotropy Haraguchi et al. (2011); Maruyama et al. (2009); Nakamura 
et al. (2009); Niranjan et al. (2010); Tsujikawa & Oda (2009). The magnetoelectric effect of 
free standing metal films has been investigated for bcc Fe(001), fcc Ni(001), and hcp Co(0001) 
subject to a uniform electric field, where the charge screening induced spin-imbalance gives 
rise to surface magnetoelectric coupling coefficients of the order of 2 — 3 x 10714 Oe cm? V71 
Duan et al. (2008), see Table 1. The magnetoelectric effect due to screening is enhanced in 
the presence of a dielectric due to the larger dielectric constant (compared to vacuum), in 
addition to effects induced by chemical bonding. For Fe/BaTiO3(001), the results of density 
functional theory (DFT) calculations predict a large surface magnetoelectric response whose 
origin is largely attributed to changes in the chemical bonding at the interface, in particular, 
to hybridization between the Ti 3d, Fe 3d and O 2p orbital states, which is found to change 
significantly upon reversal of the ferroelectric polarization direction due to the displacement 
of the Ti atoms. The magnetoelectric coupling coefficient for 2 ML Fe/BaTiO3 is estimated to 
be x = 0.01 Oe cm V~! Duan et al. (2006) while for 1 ML Fe/PbTiO3, « = 0.073 Oe cm V7! 
Fechner et al. (2008); the effect of Fe oxidation has been predicted not to affect significantly 
the magnetoelectric coupling Fechner et al. (2009). In Fe/MgO(001), a surface magnetoelectric 
coupling a; = 1.1 x 10713 Oe cm? V`! is found, significantly larger by a factor of 3.8 than 
that of a free standing Fe layer Niranjan et al. (2010); also in this system a linear change in 
the magnetocrystalline anisotropy with the applied electric field is reported. A different type 
of electric field control of magnetism consists of turning the magnetic state on and off; this 
approach has been investigated in the NiCu alloy at the composition corresponding to the 
boundary between magnetic and paramagnetic states. By capacitively charging the system, 
magnetic order can be modulated by effectively controlling the balance between the kinetic 
and exchange energies that determine the onset of magnetism Ovchinnikov & Wang (2008). 
The voltage sensitivity of ferromagnetic metallic systems near the critical temperature has also 
been studied in detail by the same authors Ovchinnikov & Wang (2009a;b). A related effect has 
been predicted for Pd, which is known to be a paramagnetic system with a Stoner parameter 
slightly short of fulfilling the condition for ferromagnetism by about 5-10%. First principles 
calculations suggest that depleting the Pd interface of charge carriers by means of an applied 
electric field can bring the Fermi level down and increase the density of states to favor an 
exchange-split (magnetic) state Kudasov & Korshunov (2007); Sun et al. (2010). A scheme 
for making the interfacial magnetoelectric effect additive consists of breaking the symmetry 
of ferromagnetic/dielectric bilayer structures by adding a non-magnetic metal at the other 
interface of the dielectric; this has been proposed by Cai et al. (2009), who have carried out 
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ab initio calculations for Fe, Co, Ni, CrO2/BaTiO3/Pt systems to demonstrate this procedure 
(see Table 1). This approach ensures that spin accumulation at one interface is not canceled 
by depletion at the other interface and that the magnetoelectric response will increase linearly 
with the number of interfaces in a superlattice structure. 


2.3 Electrostatic control of magnetism in complex oxides 

The bulk of the experimental work aiming at controlling the magnetic state electrostatically in 
complex oxide materials has focused on the “colossal” magnetoresistive (CMR) manganites, 
which are characterized by rich magnetic and electronic phase diagrams as a function of 
chemical doping (see Section 3). In these compounds, the magnetic critical temperature is 
found to coincide approximately with a peak in the resistivity versus temperature curve, 
corresponding to a metal to insulator transition Urushibara et al. (1995). In turn, the 
temperature at which the resistivity peaks has been taken as a measure of the magnetic 
ordering temperature, providing a convenient, if indirect, method of probing the magnetic 
properties of thin films and device structures. Large changes in the resistivity peak 
temperature are observed in Pb(Zro.2Tio.8)O3/Lao.8Sr0.2 MnO; (PZT /LSMO) heterostructures, 
by 35 K for a 4.0 nm LSMO film Hong et al. (2003) and 50 K for a 3.8 nm LSMO film Hong 
et al. (2005); changes up to 43 K have also been reported for 7 u.c. Lag.7Sr9.3MnO3 side-gated 
channels Pallecchi et al. (2008); and for a 5 nm Lag.gCao,2MnO3 /electrolyte field effect device, 
where a change of 32 K in the resistivity peak is observed between depletion and accumulation 
states. However, although closely related to the onset of magnetic order, the resistivity 
peak need not coincide with the magnetic critical temperature Bertacco et al. (2005); Lofland 
et al. (1997); indeed, a direct comparison between the resistivity and magnetization data in 
PZT/LSMO structures shows that while the peak in resistivity changes by 35 K between 
depletion and accumulation states, the critical temperature changes only by 20 K Vaz et al. 
(2010b). Direct measurements of the magnetic order parameter as a function of the applied 
electric field have been reported for PZT/10 nm Lag.g5Bap;5MnO3 Kanki et al. (2006), in 
PZT/4 nm Lag. g5r9.2MnO3 Molegraaf et al. (2009); Vaz et al. (2010b) structures, and in 4-6 nm 
Lag.67Sr0,33MnO3/SrTiO3 field effect device structures Brivio et al. (2010) using magnetooptic 
magnetometry. In the latter study, both top and bottom gated structures were investigated; 
the experimental results show that the critical temperature is modulated only for the top gated 
structure, i.e., for the structure where the screening occurs at the top LSMO interface Brivio 
et al. (2010). This result indicates that the properties of the bottom interface, characterized by 
the presence of an electric and magnetically “dead layer,” is less amenable to field modulation. 
For the PZT/LBMO device structures Kanki et al. (2006), a change in the magnetic hysteresis 
loops is observed as the PZT polarization is switched, with the magnetic signal decreasing 
when going from the accumulation to the depletion state, a trend opposite to that found in 
the PZT/LSMO system Molegraaf et al. (2009); Vaz et al. (2010b), whose discussion we defer 
to Section 3. 

The study of the magnetoelectric coupling in complex oxide heterostructures has also been 
carried out from the vantage point of first principles calculations. For Fe304/BaTiO3(001), 
a magnetoelectric response of a; ~ 2 x 1071? Oe cm? V~! for TiOz-terminated BaTiO3(001) 
is found Niranjan et al. (2008), comparable to the value obtained for Fe/BaTiO3 Duan et al. 
(2006); the magnetoelectric effect in both these systems is attributed to changes in the bonding 
length of the Fe cations as a function of the direction of the ferroelectric polarization, giving 
rise to large changes in the magnetic moment of Fe. A magnetoelectric coupling based 
on charge screening and in the enhancement of the spin imbalance at the Fermi level has 
been studied in a symmetric SrTiO3/SrRuOQ3/SrTiO3 structure, where the spin imbalance is 
described in terms of a spin capacitor effect, with the spin asymmetry stored at the interfaces 
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in a fashion similar to that of charge in a normal capacitor Rondinelli et al. (2008). A 
measure of the spin response of the interface is given in terms of the ratio of the surface 
spin polarization to the surface charge density, which is found to attain the value 7 = 0.37; 
this value remains the same when the SrTiO} is replaced by BaTiO3, although the change 
in magnetic moment in the SrRuO3 is much larger with BaTiO3, a direct consequence of the 
larger amount of charge required to screen the ferroelectric polarization. The magnetoelectric 
response of SrRuO3/BaTiO3(001) has also been studied, showing that, when subject to an 
electric field, the magnetic moment and the exchange splitting of SrRuO3 are modified 
due to screening, resulting in a magnetoelectric coefficient & = 2.3 x 1071? Oe cm? V-1 
Niranjan et al. (2009). For half-metallic systems, the charge carrier density also determines the 
magnetic moment, and a simple argument shows that a universal magnetoelectric constant 
due to charge screening is expected, as = pp/ ec? = 6.44 x 10714 Oe cm? V~! Cai et al. 
(2009); Duan et al. (2009). A different magnetoelectric coupling mechanism is predicted 
theoretically in Lay_,A,MnO3/BaTiO3(001) (A = Ca, Sr, or Ba) Burton & Tsymbal (2009), 
which relies on charge-induced modifications of the magnetic ground state of the CMR 
manganites Molegraaf et al. (2009). By choosing a doping level near the boundary between the 
ferromagnetic and antiferromagnetic state of Laı_,AxMnO;, the electrostatic doping arising 
from the ferroelectric polarization acts to favor either the ferromagnetic state (depletion 
state) or antiferromagnetic state (accumulation state), giving rise to a large change in the 
total magnetic moment for the two states of the ferroelectric polarization. This mechanism 
to magnetoelectric coupling had been suggested to occur in PZT/LSMO heterostructures 
Molegraaf et al. (2009) and will be discussed in more detail from an experimental perspective 
in the next section. 


3. Magnetoelectric coupling in PZT/LSMO multiferroic heterostructures 


The doped lanthanum manganites are complex oxides characterized by a strong interplay 
between charge, spin, and crystal lattice distortions, which is at the origin of the 
multifunctional behavior that is a hallmark of this class of compounds Moreo et al. (1999); 
Tokura & Tomioka (1999). Examples of the rich electronic and magnetic behavior found 
in the doped lanthanum manganites include magnetic and charge-ordered states, colossal 
magnetoresistance (CMR), high spin polarizations, and various electron transport regimes. 
The doped manganites crystallize in the pseudo-cubic AMnO; perovskite structure, where the 
12-fold coordinated A-site cations are occupied by a large ion (e.g., alkaline and rare earths), 
while the Mn cations occupy octahedrally coordinated sites Johnsson & Lemmens (2007). 
Starting with the parent compound lanthanum manganite, LaMnO;, where the Mn cations 
are in a trivalent state, the chemical substitution of La by a divalent alkaline earth removes the 
eg electron from the Mn cation, effectively adding a hole carrier to the system. The addition of 
carriers leads to profound modifications in the electronic and magnetic properties, resulting 
in complex phase diagrams as a function of chemical doping that include several electronic 
ground states Dagotto et al. (2001); Tokura (2006). This sensitivity to charge suggests that large 
susceptibilities to external electric fields can be attained when the system lies at the boundary 
separating two different ground states; by driving the system across the phase boundary using 
electrostatic doping, a change in the magnetic ground state of the system may be achieved. 
This approach has been explored in PZT /Lag gSro.2MnO3 heterostructures, where the LSMO 
system is near the boundary separating insulating and metallic ferromagnetic ground states; 
large changes in the magnetic properties are expected by using the ferroelectric field effect 
approach to modulate the charge carrier doping of the LSMO film. For optimal use of the 
ferroelectric field effect, the channel layer thickness should be comparable to the screening 
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length Ahn et al. (2003); for LSMO (x = 0.2), with a charge carrier concentration of the 
order of 10?! cm~’, the screening length has been estimated experimentally to be about 1 
u.c. Hong et al. (2005). Hence, the growth of complex oxide field effect devices requires 
precise control of the thickness down to the unit cell level. Such fine control can be achieved 
with molecular beam epitaxy Vaz et al. (2010d) or pulsed laser deposition Huijben et al. 
(2008), where the film growth can be monitored layer-by-layer in real time by using the 
oscillations in the intensity of reflection high energy electron diffraction patterns. In the 
following, we provide an overview of the recent work demonstrating a large magnetoelectric 
coupling in PZT/LSMO heterostructures, as determined by probing directly the magnetic 
order parameter using local magnetooptic Kerr effect magnetometry Molegraaf et al. (2009). 
By using advanced spectroscopy techniques, we show that the observed effect is electronic in 
origin, and that it results from a change in the valence state of the Mn cations with the change 
in the hole carrier density Vaz et al. (2010b). 

The sample structures consist of 250 nm Pb(Zro 2 Tio.s)O3 / Lao.sSto 2MnO3/SrTiO3(001), grown 
by a combination of molecular beam epitaxy for the LSMO film and off axis r.f. magnetron 
sputtering for the PZT layer. The LSMO film thickness is chosen to lie atthe transition between 
the insulating and metallic states, typically 10-12 u.c. for x = 0.2 doping. The structures are 
single crystalline, with atomically flat and sharp interfaces Vaz et al. (2010d); PZT/LSMO films 
are deposited on both unpatterned and prepatterned SrTiO3(001) substrates; the latter consist 
of Hall bar device structures, defined prior to film deposition by optical lithography and with 
dimensions optimised for optical spectroscopy measurements Vaz et al. (2010d) (see Fig. 1(b), 
inset). A 10 nm Au gate electrode is then deposited onto the PZT layer, defining the active area 
of the device (i.e, the sample region where the PZT polarization is switched), using the LSMO 
layer as the bottom contact. The Au layer is chosen to be sufficiently thin to allow transmission 
of visible light for magnetooptic Kerr effect (MOKE) magnetometry measurements. MOKE 
relies on the fact that the polarization of light is modified upon reflection from a magnetic 
surface; it is a technique particularly well suited for this study, since it allows a direct and 
local measurement of the magnetic order parameter. In one implementation of this technique, 
a linearly polarized laser beam is reflected off the sample surface (with the plane of incidence 
parallel to the applied field direction, called the longitudinal MOKE geometry), and the 
Kerr rotation or ellipticity, which is proportional to the magnetization, is measured using a 
polarimeter unit Vaz et al. (2010e). 

The individual electric and magnetic characteristics of the PZT/LSMO heterostructure are 
shown in Fig. 1(a) and (b), respectively. The electric polarization versus electric field (P-E) 
response shows abrupt electric switching and a saturation polarization of about 85 uC cm~?; 
the magnetic hysteresis (M-H) curves, for both the accumulation and depletion states, show 
that the system is ferromagnetic at 100 K and that there is a marked difference in the 
magnetic properties for the two states of the ferroelectric polarization, namely, a larger 
coercivity and a smaller saturation magnetization for the accumulation state as compared 
to the depletion state. These individual ferroic curves are the classical hysteresis curves 
of ferroelectrics and ferromagnets; the magnetoelectric coupling is demonstrated by the 
magnetic response of the system as a function of the applied electric field (M-E loop) shown in 
Fig. 1(c), where the saturation magnetization is found to switch hysteretically and reversibly 
between a low and high magnetic moment at the electric field values corresponding to 
the switching of the ferroelectric polarization. This result demonstrates the presence of a 
magnetoelectric coupling in these multiferroic heterostructures, showing in particular that the 
direction of the PZT ferroelectric polarization modifies the magnetic state of the LSMO layer. 
Note that the difference in magnetic moment persists at zero applied electric field, which 
excludes electrostrictive or piezoelectric effects (strain) as being the cause of the observed 
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magnetoelectric effect. The magnetoelectric effect in the PZT/LSMO system is robust and 
has been observed in structures where the LSMO film has been grown by off axis magnetron 
sputtering Molegraaf et al. (2009) and molecular beam epitaxy Vaz et al. (2010b). 
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Fig. 1. Composite multiferroic heterostructure characterized by a coupling between the 
classical ferroelectric (a) and magnetic (b) ferroic responses, as shown in (c). While the M-H 
and P-E loops are standard, the M-E characteristic is new. The example shown is for a 
PZT/12 u.c. Lag gSrg.2MnO3 composite multiferroic heterostructure, from Vaz et al. (2010b; 
2011). Figures (a) and (b): Copyright 2010 by The American Physical Society. Figure (c) 
reprinted with permission from C.A.F. Vaz et at., J. Appl. Phys., 109, 07D905 (2011). 
Copyright 2011, American Institute of Physics. 


The variation of the magnetization as a function of temperature for the two states of the PZT 
polarization is shown in Fig. 2(a) Molegraaf et al. (2009); Vaz et al. (2010c). The magnetization 
curves are characteristic of a ferromagnetic system, with a critical temperature separating 
the high temperature paramagnetic regime and a ferromagnetically ordered state at low 
temperatures. What is striking in these data is that the direction of the PZT polarization 
determines the magnetic properties of the system, including an increase in the critical 
temperature (by about 20 K) and a decrease in the ground state magnetization when switching 
from the depletion to the accumulation state (in agreement with the M-H hysteresis curve of 
Fig. 1(b)). One sees that the state corresponding to the highest Curie temperature has the 
lowest saturation magnetization, which agrees qualitatively with what is expected from the 
behavior of bulk LSMO, since hole doping changes the ionic state of Mn+, with spin S = 2, 
to Mn**, with spin S = 3/2, so that we expect a decrease in magnetization with increasing 
doping Jonker & van Santen (1950). 

As shown in Fig. 1(c), the magnetoelectric response of the PZT/LSMO multiferroic 
heterostructure is strongly non-linear, and the response of the system is best described 
in terms of an effective magnetoelectric susceptibility AM/AE = AM/2E., where E. is 
the ferroelectric coercive field; at 100 K, we find AM/AE = 6.2 x 10-2 Oe cm v1, 
This value is significantly larger, by 2-3 orders of magnitude, than typical magnetoelectric 
coupling coefficients of single-phase multiferroics and comparable to the value obtained for 
strain-mediated composites Fiebig (2005); Ma et al. (2011); Vaz et al. (2010a). Given the 
interfacial nature of the magnetoelectric effect in this system, one alternate measure of the 
magnetoelectric effect is given in terms of the surface (interface) magnetoelectric coefficient 
as, which is obtained by multiplying AM/2E, by the LSMO film thickness, yielding a; = 
2.9 x 107? Oe cm? V~! at 100 K. In Vaz et al. (2010c), the variation of the magnetoelectric 
response of PZT/LSMO as a function of the temperature is studied, as shown in Fig. 2(b). 
AM/2E. is found to have a strong, non-monotonic, temperature variation, including a change 
in signal at around 150 K; in particular, one finds that the magnetoelectric response is largest 
at around 180 K, with AM/AE = —13.5 x 107°? Oe cm V~! (a; = —6.3 x 107? Oe cm? V~}), 
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Fig. 2. (a) Magnetic response of the system as a function of temperature for the two states of 
the PZT polarization. (b) Variation in the magnetoelectric coupling coefficient AM/2E. as a 
function of temperature. After Vaz et al. (2010c). 


near the magnetic critical point, where the difference between the magnetization for the 
depletion and accumulation states is largest, as shown in Fig. 2(a). The change in sign of the 
magnetoelectric response can be understood with reference to the magnetization data shown 
in Fig. 2(a). One finds that 150 K separates a low temperature regime, where the magnetization 
for the depletion state is larger than that of the accumulation state, from a high temperature 
regime (up to the paramagnetic state), where the opposite behavior occurs due to the faster 
decay of the magnetization of the system in the depletion state on approaching its magnetic 
critical temperature. What is particularly promising in this system is that AM/AE peaks at 
temperatures near the critical region, showing that the largest magnetoelectric response is 
achieved at the highest temperature where the system remains magnetic. A key aspect of 
composite multiferroics is that they need not obey reciprocity relations that limit the size of 
the magnetoelectric susceptibility in single phase compounds Brown et al. (1968). 

The microscopic origin of the magnetoelectric effect is investigated by x-ray absorption near 
edge spectroscopy (XANES). XANES is a technique particularly well suited to study the 
changes in the electronic structure of Mn as a function of the applied electric field, due to 
its sensitivity to changes in the valence state. In XANES, one excites electrons from a core 
shell to empty states in the valence band (1s to 4p states for the case of the Mn K edge). 
This excitation energy depends on the formal valence state of the cation, with the more 
electronegative cations requiring higher energy photons to excite core electrons Kirichok et al. 
(1985). One key advantage of x-ray absorption spectroscopy is its ability to probe buried 
layers, which is not possible with some other spectroscopic techniques, such as photoelectron 
spectroscopy. Details of the experimental set-up and measurement conditions are given in 
Vaz et al. (2010b;e). 

The room temperature variation of the x-ray light absorption as a function of the incident 
photon energy across the absorption K edge of Mn, for the depletion and accumulation states, 
is shown Fig. 3(a) Vaz et al. (2010b). The key result is the observation of an energy shift in 
the absorption edge of Mn between the depletion and accumulation states, by 0.3 eV, which 
corresponds to about 10% of the total shift expected for the full LSMO doping range Bindu 
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Fig. 3. (a) Room temperature XANES results for the two polarization states of the PZT. (b) 
Difference in x-ray absorption for the two PZT polarization states; the full line models this 
difference assuming a rigid shift in the Mn absorption edge. (c) Variation of the x-ray light 
absorption as a function of the applied gate voltage at a fixed energy, E = 6549.7 eV. The 
error bars reflect counting statistics. After Vaz et al. (2010b). 


et al. (2005); Shibata et al. (2003). The shift in the absorption edge can also be observed 
in the difference x-ray absorption spectrum, shown in Fig. 3(b). The observation of a shift 
in the absorption K edge of Mn shows that the valence state of Mn changes with charge 
doping induced by the ferroelectric polarization and demonstrates the electronic nature of the 
magnetoelectric coupling in these multiferroic heterostructures. The same result is observed 
at low temperature (20 K). From the observed energy shift we estimate the average change 
in valency across the LSMO film, Ax = 0.1 per Mn, using the results by Shibata et al. 
(2003) showing a linear variation in the energy shift with the formal average valency of Mn, 
AE = 3.0x, where x is the LSMO doping. 

One can now compare the change in valency as determined by XANES with the values 
obtained from the electric and magnetic characterization of the same device structure. From 
the saturation electric polarization, P; = 85 uC cm~?, one can deduce the expected change 
in the carrier doping, An = 0.13 e/Mn. This value is in good agreement with the change in 
valency determined from XANES, Ax = 0.1 per Mn, showing that the electric polarization is 
screened effectively by hole carriers from the LSMO layer. Changing the valency of Mn from 
3+ to 4+ leads to a change in spin state by 1/2 and to a change in magnetic moment of gS ~ 1 
Up (where g ~ 2 is the g-factor), so one expects a change in the average magnetic moment 
of 0.1 ug/Mn. However, from MOKE (Fig. 1(b)) in combination with SQUID magnetometry, 
one obtains a change in magnetic moment of about 0.76 up/Mn. This change is much larger 
than the expected change of 0.1 #g/Mn from the change in spin state. One sees therefore that 
the change in magnetic moment cannot be explained simply by the change in the Mn spin 
state. From these results, and from the amplitude of the change in the magnetic moment, one 
concludes that a change in both the spin state and spin configuration must occur, whereby 
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the depletion state remains in a ferromagnetic state, while the accumulation state exhibits an 
antiferromagnetic configuration confined to the PZT interface where charge screening occurs. 
A model depicting the magnetic reordering is shown in Fig. 4. In this picture, the interface 
layer in the accumulation state consists of strongly depopulated 3d z? states, weakening 
the double-exchange interaction at these orbitals and favoring superexchange, leading to an 
antiferromagnetic coupling of the interfacial layer. This picture is in qualitative agreement 
with first principles calculations in strained LSMO Fang et al. (2000), and in particular with 
recent ab initio calculations by Burton and Tsymbal on La;_,BayMnO3/BaTiO3 Burton & 
Tsymbal (2009), where the low energy state of the system at x = 0.5 in the accumulation state 
corresponds to an antiparallel spin configuration of the interfacial spins. In this simplified 
model, the change in the interfacial spin configuration leads to a change in magnetic moment 
of about 0.6 ug/Mn between depletion and accumulation states, which agrees with the 
experimental results. This mechanism gives rise to a much more dramatic change in the 
average magnetic moment and explains the large magnetoelectric coupling observed in this 
system. One direct consequence of these results is that one can control electrostatically 
both the spin state and magnetic spin configuration of the interfacial LSMO layer. Given 
that this layer is also responsible for charge transport, and given the high spin polarization 
characteristic of LSMO Park et al. (1998), this system constitutes a potential platform for the 
design of spintronics devices, such as lateral spin-valve structures, where the magnetic and 
transport properties are controlled by applied electric fields. 


PZT 


depletion state accumulation state 


Fig. 4. Schematic model of the spin configurations in LSMO at the PZT interface for the 
depletion and accumulation states, showing the changes in the Mn and O orbital states and 
the expected changes in the magnetic moment per layer. The arrows indicate the spin 
orientation in the Mn cations and n denotes the unit cell number below the PZT. The Mn d 
orbitals are drawn in orange, and gray, and the lobes of the p orbitals are shown around the 
oxygen atoms (red). After Vaz et al. (2010b), Copyright 2010 by The American Physical 
Society. 


4. Conclusion and outlook 


To summarize, we illustrate how novel functionalities can be engineered by exploring the 
new phenomena that arise at the interface between dissimilar materials and how such 
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an approach can be used to achieve electrostatic control of magnetism in multiferroic 
heterostructures. By exploiting the sensitivity of the ground state properties of the CMR 
manganites to charge, we showed that large, charge-driven, magnetoelectric coupling in 
PZT/LSMO multiferroic heterostructures can be achieved. The effect is electronic in origin, as 
demonstrated by advanced spectroscopic techniques, and is therefore compatible with current 
CMOS technology, in particular since the growth of crystalline oxides on Si(001) is now well 
established McKee et al. (1998; 2001); Reiner et al. (2008). Further, we show that both the 
spin state and the magnetic configuration can be controlled electrostatically. One research 
challenge remains the optimization of the materials properties to allow room temperature 
operation; possible solutions are the use of optimal doping or the use of other magnetic oxides, 
such as the double perovskites. Another area that needs to be addressed is the dynamic 
magnetoelectric response of these multiferroic heterostructures, both in terms of the electric 
field modulation of the spin wave spectra and in terms of the electrostatic switching time of 
the magnetic state. Given that in these multiferroics the processes are electronically driven, 
one may expect very fast switching times, limited by the ferroelectric switching process. As 
this chapter aims to show, demonstrating the electrostatic switching of magnetism in a range 
of materials has been achieved as a proof of concept, and much remains to be investigated 
in terms of materials optimization, development and characterization of new systems, and 
understanding of the interfacial processes mediating the magnetoelectric coupling. The work 
achieved thus far illustrates how novel device structures enabling electric field control of the 
magnetic spin state can be tailored to provide the requisite characteristics that will make them 
potential candidates for next generation electronic devices. 
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